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Abstract. The ROSAT Deep Survey in the Lockman 
Hole is the most sensitive X-ray survey performed to 
date, encompassing an exposure time of 207 ksec with 
the PSPC and a total of 1.32 Msec with the HRI aboard 
ROSAT. Here we present the complete catalogue of 50 
X-ray sources with PSPC fluxes (0.5-2 keV) above 5.5 x 
10 -15 erg cm~ 2 s _1 . The optical identifications are dis- 
cussed in an accompanying paper (Schmidt et al., 1997). 
We also derive a new log(N)-log(S) function reaching a 
source density of 970 ± 150 deg~ 2 at a limiting flux of 



10 



-15 



erg cm 



At this level 70-80% of the 0.5-2 keV 



X-ray background is resolved into discrete sources. Utiliz- 
ing extensive simulations of artificial PSPC and HRI fields 
we discuss in detail the effects of source confusion and in- 
completeness both on source counts and on optical identi- 
fications. Based on these simulations we set conservative 
limits on flux and on off- axis angles, which guarantee a 
high reliability of the catalogue. We also present simula- 
tions of shallower fields and show that surveys, which are 
based on PSPC exposures longer than 50 ksec, become 
severely confusion limited typically a factor of 2 above 
their 4cr detection threshold. This has consequences for 
recent claims of a possible new source population emerg- 
ing at the faintest X-ray fluxes. 
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1. Introduction 

The study of the X-ray background (XRB), originally dis- 
covered 35 years ago (Giacconi et al., 1962), with imaging 
X-ray telescopes has progressed rapidly in the last few 
years. After a pioneering start with deep Einstein expo- 
sures (Giacconi et al., 1979, Griffiths et al., 1983, Hamil- 
ton et al., 1991), ROSAT deep survey observations were 
able to resolve the majority of the soft XRB into discrete 
sources (Hasinger et al., 1993, hereafter H93). First deep 
imaging exposures with ASCA have extended these stud- 
ies into the hard X-ray band (Inoue et al. 1996, Cagnoni 
et al., 1997). The X-ray band is one of the few regions of 
the electromagnetic spectrum, where the integrated emis- 
sion of discrete sources, summed over all cosmic epochs, 
is dominating the celestial emission. The XRB therefore 
provides a strong constraint for cosmological evolution 
models (see e.g. Comastri et al., 1995; Zdziarski 1996). 
In order to understand the populations contributing to 
the X-ray background and their cosmological evolution it 
is necessary to obtain complete optical identification of 
the sources detected in X-ray deep surveys. The major- 
ity of optically identified X-ray sources at faint fluxes are 
active galactic nuclei (AGN), i.e. QSOs and Seyfert galax- 
ies with broad emission lines (Shanks et al., 1991; Boyle 
et al., 1993; Bower et al, 1996; Page et al., 1996). Using 
the large samples of X-ray selected AGN now becoming 
available, their luminosity function and its evolution with 
redshift can be studied in detail (Boyle et al., 1995; Page 
et al., 1996; Jones et al., 1996). The current models, which 
are based on pure luminosity evolution, predict that the 
broad-line AGN contribute only 30-50% to the soft X-ray 
background and much less to the hard X-ray background. 
Postulating a large population of faint, intrinsically ab- 
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sorbed AGN, some models predict a much larger fraction 
of the soft XRB, and possibly all of the hard XRB as due 
to AGN (see e.g. Comastri et al., 1995; Zdziarski 1996). 

There are, however, reports that at the faintest X-ray 
fluxes a new population of X-ray active, but optically rela- 
tively normal narrow-emission line galaxies (NELG) start 
to dominate the X-ray counts, and could ultimately con- 
tribute the majority of the XRB (Georgantopoulos et al. 
1996, Griffiths et al. 1996, McHardy et al., 1997). These 
claims are based on deep ROSAT PSPC surveys which - 
as we will show - start to be severely confusion limited at 
the faintest fluxes. Care has therefore to be taken to assess 
the systematic position errors and misidentifiction likeli- 
hood at faint fluxes. In our own deep survey identification 
work we have therefore taken measures to minimize confu- 
sion and to increase the reliability of the X-ray catalogue. 
We restricted the sample to X-ray fluxes well above the 
detection threshold at a reliable flux level derived from de- 
tailed simulations. Secondly, we invested massive amounts 
of ROSAT HRI time to cover our deepest survey field in 
the Lockman Hole, this way minimizing position errors 
and confusion. In this paper (paper I) we describe the 
derivation of the X-ray source catalogue in the Lockman 
Hole. Section 2 gives a summary of the observations, Sect. 
3 describes the new maximum likelihood (ML) crowded- 
field detection algorithm developed for deep X-ray survey 
work, Sect. 4 discusses the necessary astrometric correc- 
tions for ROSAT PSPC and HRI images and Sect. 5 gives 
the final source catalogue. The verification of the analysis 
procedure through simulations and the derivation of a fi- 
nal log(N)-log(S) function are discussed in Sect. 6 and a 
discussion of the results is given in Sect. 7. The complete 
optical identification of the X-ray source catalogue is de- 
scribed in paper II (Schmidt et al., 1997) and deep VLA 
radio observations of the field are discussed by deRuiter 
et al. (1997). 

2. Observations 

The data presented here consists of all observations of the 
ROSAT Deep Survey accumulated in the period 1990- 
1997 in the direction of the Lockman Hole, which is one 
of the areas on the sky with a minimum of the galactic 
Hydrogen column density; Nh is roughly 5 x 10 19 cmT 2 
in this field (Lockman ct al., 1986). Table 1 gives a sum- 
mary of these observations. Specifically, a total exposure 
time of 207.41 ksec was accumulated with PSPC (Pfeffer- 
mann and Briel 1992) pointings centered at the direction 
i?A(2000) = 10' l 52 m , L>.EC(2000) = 57°21'36". These 
data have been partially presented in H93. For the PSPC 
data all good time intervals selected by the ROSAT stan- 
dard analysis (SASS; Voges 1992) have been analysed. The 
individual datasets obtained in different semesters, which 
are affected by residual erratic aspect errors of order 10 
arcsec, were all shifted to a common celestial reference 
system using preliminary optical identifications. A PSPC 



image with a pixel size of 5 arcsec, covering the most sen- 
sitive area of the field of view (FOV), with a radius of 
~ 20 arcmin, has been accumulated in the standard PSPC 
"hard" energy band 0.5-2 keV (pulseheight channel 52- 
201). This choice optimizes the angular resolution of the 
PSPC and the signal to noise ratio for the detection of 
faint extragalactic sources. It is also relatively insensitive 
to neutral hydrogen absorption column densities below 
~ 10 21 cm~ 2 . For the calculation of hardness ratios we 
additionally analysed the PSPC data in the "soft" band 
and further subdivided the hard band into two: HI and 
H2. Details about the energy bands are given in Table 2. 
Fig. |l|a shows a contour plot of the PSPC image in the 
hard band. 

The same area has also been covered by a set of raster 
scan observations with the ROSAT HRI (David et al., 
1996). A total of ~ 100 pointings of roughly 2 ksec ex- 
posure each, has been placed on a regular grid with a 
step size of ~ 2 arcmin. The HRI FOV has a size of 
roughly 36 x 36 arcmin. The inner, most sensitive part 
of the field, where the off-axis blur of the telescope can 
be neglected, has a radius of roughly 8 arcmin. The raster 
scan step size is much smaller than the HRI FOV, thus 
the resulting exposure covers the survey area smoothly, 
with a rather homogeneous point spread function (PSF) 
across the field. The total exposure time accumulated for 
the raster scan is 205.50 ksec, comparable to the PSPC 
pointing. The individual HRI pointing data provided by 
the ROSAT Standard Analysis System SASS (Voges 1992) 
have been merged into a single mosaique, after separate 
astrometrical solutions have been applied to each dataset, 
in order to correct for residual boresight and scale fac- 
tor errors in the HRI data (see below). The HRI raster 
image has been accumulated in the restricted pulseheight 
channel range 1-9, which significantly reduces the detec- 
tor background with a minimum loss (~ 7%) of cosmic 
X-ray photons (David et al., 1996). Using this channel se- 
lection, the HRI is sensitive to cosmic X-ray photons in 
the full ROSAT energy range 0.1-2.4 keV. Fig. |l|b shows 
the resulting HRI raster scan image on the same scale as 
the PSPC data. 

From Fig. [I] it is obvious, that the HRI raster scan, al- 
though having significantly higher resolution, is not nearly 
as sensitive as the PSPC observation which, however, is 
reaching the confusion limit (see below) . In an attempt to 
reach an unconfused sensitivity limit comparable or deeper 
than the PSPC exposure, at least in a substantial fraction 
of the PSPC FOV, we performed an ultradeep HRI expo- 
sure of 1 Msec in a single pointing direction. The ROSAT 
attitude system shows residual pointing errors on the or- 
der of several arcsec, which can lead to a corresponding 
positional shift of the resulting image (see below). Uncor- 
rected aspect errors can also lead to a substantial ellip- 
ticity of the HRI images of point sources (David et al., 
1996), thus reducing resolution and sensitivity. To some 
degree these aspect errors can be corrected for, if one or 
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Table 1. Observation Summary 



Start Date 


End Date 


Instrument 


Exp. [s] 


1990 Apr 16 


1991 May 21 


PSPC 


67989 


1991 Apr 25 


1991 May 5 


HRI raster 


86778 


1991 Oct 25 


1991 Nov 2 


PSPC 


24327 


1991 Oct 27 


1991 Nov 9 


HRI raster 


116635 


1992 Apr 15 


1992 Apr 24 


PSPC 


66272 


1992 Apr 18 


1992 Apr 18 


HRI raster 


2082 


1992 Nov 29 


1992 Nov 29 


PSPC 


2082 


1993 Apr 26 


1993 May 9 


PSPC 


46740 


1994 Oct 21 


1994 Nov 7 


HRI 


106565 


1995 Apr 15 


1995 May 11 


HRI 


294483 


1995 Oct 26 


1995 Nov 27 


HRI 


200358 


1996 May 1 


1996 May 29 


HRI 


319590 


1997 Apr 15 


1997 Apr 28 


HRI 


191094 



more bright optically identified X-ray sources are present 
in the image. In order to allow this procedure, we selected 
a pointing direction for the ultradeep HRI exposure, which 
is inside the PSPC field of view, but shifted about 10 
arcmin to the North-East of the PSPC center, this way 
covering a region containing about 10 relatively bright X— 
ray sources known from the PSPC and the HRI raster 
scan. The ultradeep HRI exposure is centered at the direc- 
tion iL4(2000) = 10 h 52 m 43 s , DEC(2000) = 57°28'48". In 
some situations the automatic ROSAT HRI standard anal- 
ysis system is too conservative in rejecting time intervals 
with supposedly high background. In particular using the 
background-optimized HRI channel range 1-9, we found 
an optimum point-source detection sensitivity including 
all exposure intervals with good attitude data. For the ul- 
tradeep survey we therefore selected events using a custom 
made procedure after having joined the SASS accepted 
and rejected event files. This resulted in a net observing 
time of 1112.09 ksec, a gain of ~ 20% compared to the 
SASS products. The individual datasets obtained in differ- 
ent semesters were shifted to the same reference system as 
the PSPC data, using preliminary optical identifications. 
The HRI image in the pulseheight channel range 1-9 is 
shown in Fig. |l|c on the same scale as the PSPC image, 
which makes the shift to the North-East quite apparent. 

3. Detection algorithm 

The images in Fig. [j] have been analysed using an im- 
proved version of the interactive analysis system EXSAS 
(Zimmermann et al. 1994). The first steps of the detec- 
tion procedure, i.e. the local detect algorithm LDETECT, 
the bicubic spline fit to the background map and the map 
detection algorithm MDETECT have been described in 
H93. The detection threshold of these algorithms has been 
set at a very low likelihood value so that the resulting 
list of possible source positions contains several hundred 
spurious detections. This position list has then been fed 
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Fig. 1. ROSAT Deep Survey images in the Lockman Hole: (a) 
PSPC 207 ksec pointing, (b) HRI 205 ksec raster scan and (c) 
HRI 1112 ksec pointing. The images are all centered at the 
PSPC pointing direction (see text). 
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into the "multi-ML" crowded field detection and param- 
eter estimation algorithm first described by Hasinger et 
al. (1994a). This method, which has been implemented 
into EXSAS, works on binned image data and fits su- 
perpositions of several PSF profiles on top of an external 
background to sections of the data (typically sub-images 
of arcmin size). Best-fit positions and normalizations are 
obtained by maximizing the likelihood statistic C (Cash 
1979) or, correspondingly, by minimizing the quantity 

S = -2lnC = -2j2(Ymod(i,j) - N(i,j)lnY mod (i,j) (1) 



^rnod ( 

i, j) is the sum of all model point source contribu- 
tions plus the background value in the image pixel [i,j] and 
N is the measured number of photons in pixel [i,j]. For 
ROSAT PSPC and HRI data the multi-component PSF 
(Hasinger et al., 1994b; David et al., 1996) is approximated 
by a single Gaussian function with a width increasing with 
off-axis angle. The significance AS of any of the individ- 
ual point source contributions is then estimated by a like- 
lihood ratio test between the best fit with and without the 
corresponding component. If AS falls below a threshold 
value, the corresponding component is omitted from the 
next iteration and only significant components are main- 
tained. Since AS follows a \ 2 distribution, the errors of 
the best-fit parameters (68% confidence single parameter 
of interest) are determined by searching the parameter 
space for which AS = S — S m i n = 1. Before the next re- 
gion of interest is fit, the current PSF contributions are 
added to the background model. 

The complete source detection and parameter esti- 
mation procedure is too complicated to validate its re- 
sults and e.g. quantify its detection efficiency analytically. 
We therefore performed extensive simulations of artificial 
ROSAT fields run through the same set of algorithms to 
study their properties (see below). Based on these sim- 
ulations we have chosen a likelihood detection threshold 
of AS = 16, corresponding to a 4ct detection. Formally, 
the probability of a spurious detection with such a thresh- 
old is 6 x 10~ 5 per resolution element, which would cor- 
respond to less than one spurious detection per field of 
view. We will see in Sect. 6 that in practice the effective 
detection limit is set by source confusion. Observed fluxes 
in the 0.5-2 keV band are calculated for each source us- 
ing the detected counts, the corresponding exposure time 
and vignetting correction, an energy-to-flux conversion 
factor (ECF) and a PSF-loss correction factor (PCF) de- 
termined from simulations (see below). In order to derive 
the ECF we assumed a power law spectrum with photon 
index 2 and galactic absorption, folded through the instru- 
ment response. Table 2 gives a summary of the various en- 
ergy bands and correction factors for different detectors. 
For the PSPC hard band, because of the restricted energy 
range, the ECF is almost insensitive to Nh and to the as- 
sumed spectral shape for a wide range of absorptions and 



spectral indices. On the contrary, the HRI is very sensi- 
tive to absorption and spectral shape because of its lack of 
energy resolution. Unabsorbed sources with steep spectra 
are favoured by the HRI. 



Table 2. Energy bands and correction factors 



Detector 


Band 


Energy 


Pulseheight 


ECF" 


PCF" 






[keV] 


Channels 






PSPC 


H 


0.5-2.0 


52-201 


0.836 


0.90 


PSPC 


S 


0.1-0.4 


11-41 


1.519 


0.90 


PSPC 


HI 


0.5-0.9 


52-90 


0.349 


0.90 


PSPC 


H2 


0.9-2.0 


91-201 


0.487 


0.90 


HRI 




0.1-2.4 


1-9 


0.586 


0.92 



a Energy to counts conversion factor in cts/s for a source 
with 0.5-2 keV flux 10~ n erg cm' 2 s' 1 

b PSF loss correction factor from simulations 



For the calculation of hardness ratios, PSPC images 
were also analysed in the soft (S) and the split hard (HI, 
H2) bands (see Table 2), however with fixed positions. The 
standard ROSAT soft and hard hardness ratios HRI and 
HR2 have been used. 



HRI = 



H-S 
H + S 



HR2 = 



H2-H1 



H2 + H1 



(2) 



4. Astrometric corrections 



The ROSAT SASS performs various corrections to the in- 
dividual PSPC and HRI photons, including detector lin- 
earisation, optical distortion, borcsight and attitude cor- 
rection. The large set of X-ray sources with reliable opti- 
cal identifications in the Lockman Hole allows for the first 
time to check and crosscalibrate with high accuracy the 
ROSAT astrometry with absolute celestial positions. 

A complete sample of X-ray sources in the Lockman 
Hole with fluxes down to 5.5 x 1CP 15 erg cm' 2 s' 1 has 
been optically identified (Schmidt et al., 1997, paper II). 
The positions of the optical counterparts have been de- 
rived from Palomar 200" CCD drift scans and have an 
absolute accuracy of <~ 0.5 arcsec. Using the set of point- 
like X-ray sources optically identified with AGN and stars, 
we could derive astrometrical solutions for the PSPC and 
HRI data, respectively. For the deep PSPC and HRI point- 
ing observations astrometrical solutions were fit to the de- 
viation between X-ray and optical positions. The trans- 
formations contained the following parameters: A(RA) 
and A(DEC) position shifts, A($) roll angle and F(Sc) 
scale factor correction. The best-fit parameters were de- 
termined by a \ 2 ~^i using the individual statistical po- 
sition errors of the X-ray sources and assuming an error 
of 0.5" for the optical positions. In order to take account 
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5 



Instrument 


Nr 


J-J sys 


X 


V 


c ) 


A 

L\ 




) 


At R A\ 


At n pr 1 ! 

L\ I U Hi 1 




IDs 


["] 










I'] 




["] 


["] 


PSPC Pointing 


47 


3.0 


1.28 


1.0020 ± 


0.0007 


-2.9 


± 


2.0 


-0.2±0.5 


-0.1±0.5 


HRI Pointing 


32 


0.5 


1.05 


0.9972 ± 


0.0006 


1.2 


± 


1.4 


-0.1±0.3 


0.3±0.3 


HRI Raster 


31 


0.5 


1.07 


1.0004 ± 


0.0006 


1.0 


± 


1.5 


-0.3±0.4 


-0.3±0.4 


Combined 


43 


0.0 


1.36 


1.0001 ± 


0.0003 


-1.1 


± 


0.7 


-0.1±0.2 


0.0±0.2 



of possible additional deviations in the X-ray positions, 
which may result from residuals in the detector lineariza- 
tion or from X-ray source confusion, we added in quadra- 
ture a constant systematic position error E sys to the sta- 
tistical error for each source. The size of this additional 
error (see column 3 of Table 3) was chosen such that the 
reduced x 2 approached a value of 1 . 

Table 3 shows the result of the astrometric solutions. 
As expected, the residual position shifts are all negligi- 
ble, because all data have been shifted to a common sys- 
tem before. The roll angle correction applied by the SASS 
is already very good; the roll angle corrections A(<I>) in 
our deep datasets are all smaller than a few arcmin. The 
only significant astrometrical correction which is neces- 
sary for both PSPC and HRI is a scale factor correction 
of order 0.2-0.3%. For the PSPC this correction factor is 
marginally significant and can probably be neglected for 
most datasets. For the HRI the scale factor error is con- 
sistent with the estimate derived from HRI data of M31 
(David et al., 1996), but our data is more accurate and 
provide the first significant measurement of this quantity. 
The HRI scale factor error corresponds to a position de- 
viation of about 4 arcsec across the HRI FOV and needs 
to be taken into account to achieve highest astrometrical 
accuracy. The residual systematic position errors across 
the HRI FOV are very small (0.5"), while the PSPC has 
residual systematic errors of order 3" , most likely because 
of source confusion and the broader PSF (see below). 

For the HRI raster scan data, every individual short 
pointing was first shifted to a common celestial reference 
system and corrected for the above scale factor error, be- 
fore all photons were merged to one dataset. As a test 
we applied the same astrometrical transformation fit as to 
the pointed observations. The fact that all fit parameters 
are consistent with their expected values and that the re- 
quired residual systematic errors are of the same order as 
for the individual HRI pointing shows that the shift- and 
merge procedure was successful. Finally, the last row in 
Table 3 shows an astrometric fit to the X-ray source cat- 
alogue combined from all three datasets (see below). As 
expected, all astrometric parameters are consistent with 
their expected values. The reduced x 2 of the combined 
fit is acceptable without the addition of extra systematic 
errors. The distribution of positional errors for uniquely 



identified sources is therefore as expected, thus validating 
the calculation of the positions and errors. 

5. X ray source catalogue 

The scale factor corrections determined in the previous 
section were applied to the three individual detection lists 
from PSPC pointing, HRI pointing and HRI raster scan 
observation and the corresponding systematic position er- 
rors were added in quadrature to the statistical errors de- 
termined in the multi-ML detection procedure. Finally, 
the three individual position lists were merged into a sin- 
gle source list. For all detections, which were positionally 
coincident within their 90% error radii, weighted average 
positions and position errors were calculated. These form 
the basis for the final X-ray source catalogue, presented 
in Table 4. This catalogue contains all objects with PSPC 
fluxes higher than 1.11 x 10~ 14 erg cmT 2 s^ 1 at off-axis 
angles smaller than 18.5 arcmin and all objects with fluxes 
larger than 0.56 x 10~ 14 erg cm~ 2 s _1 inside an off-axis 
angle of 12.5 arcmin. 

Column (1) of Table 4 gives the official ROSAT name, 
column (2) an internal source number. Columns (3) and 
(4) give the weighted average coordinates of the X-ray 
sources for an equinox of J2000.0. Column (5) shows the 
la position error of the source, including statistical and 
systematic errors. A capital P, H, or R after the position 
error indicates, whether the dominant weight to the X-ray 
position comes from the PSPC pointing, HRI pointing or 
HRI raster scan, respectively. Column 6 gives the 0.5-2 
keV flux of the source in units of 10~ 14 erg cmT 2 s~\ 
derived from the ROSAT PSPC hard band and its \a 
error. Column (7) and (8) give the soft and hard hardness 
ratios defined above. 

6. Verification of the analysis procedure by simu- 
lations 

The ROSAT deep survey exposures are probing the lim- 
its of observational and data analysis procedures. In or- 
der to obtain a reliable and quantitative characteriza- 
tion and calibration of the source detection procedure, 
detailed simulations of large numbers of artificial fields, 
analysed through exactly the same detection and param- 
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Table 4. Source Catalogue 



Source Name Nr. R.A. (2000) Dec. (2000) Err. fx HR1 HR2 



RX 


J105421. 


1+572545 


28 


10 


54 


21.1 


57 


25 


44.5 


1 


OH 


19. 


.90 


± 





,41 


1. 


.00 


± 





.06 


-0 


.27 


± 





.02 


RX 


J105131 


.1+573440 


8 


10 


51 


31.1 


57 


34 


40.4 





8H 


11. 


84 


± 





33 


-0 


.24 


± 





.02 


-0 


.47 


± 





.02 


RX 


J105316. 


.8+573552 


6 


10 


53 


16.8 


57 


35 


52.4 





9H 


8 


,82 


± 





,33 


-0 


.29 


± 





.02 


-0 


.40 


± 





.02 


RX 


J105239. 


.7+572432 


32 


10 


52 


39.7 


57 


24 


31.7 





8H 


6 


.50 


± 





,30 


-0 


.55 


± 





.02 


-0 


.16 


± 





.03 


RX 


J105335. 


.1+572542 


29 


10 


53 


35.1 


57 


25 


42.4 





8H 


5 


.13 


± 





,25 


-0. 


.34 


± 





.03 


-0 


36 


± 





.03 


RX 


J105331. 


.8+572454 


31 


10 


53 


31.8 


57 


24 


53.9 





9H 


3 


.57 


± 





18 


-0 


.18 


± 





.04 


-0 


.40 


± 





.03 


RX 


J105339. 


.7+573105 


16 


10 


53 


39.7 


57 


31 


5.3 





9H 


3 


.56 


± 


0. 


.18 


-0. 


.59 


± 





.02 


-0 


.43 


± 





.03 


RX 


J105020 


,2+571423 


56 


10 


50 


20.2 


57 


14 


22.8 


1 


.8R 


3 


.52 


± 


0. 


.17 


-0 


.44 


± 





.02 


-0 


.37 


± 





.04 


RX 


J105201. 


,5+571044 


62 


10 


52 


1.5 


57 


10 


44.2 


1 
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Fig. 2. Comparison of detected log(N)-log(S) relation to input source counts in the 0.5-2 keV band. The dotted line shows the 
differential source counts n(S) input to the simulations, normalized to a Euclidean slope (see equation 3 and 4). The filled and 
open circles show the detected differential log(N)-log(S), derived from (a) 66 simulated PSPC fields of 200 ksec exposure, and 
(b) 100 simulated HRI fields of 1 Msec exposure each. Filled symbols correspond to sources at off-axis angles smaller than 12.5 
arcmin for the PSPC and 10 arcmin for the HRI, respectively. Open symbols refer to larger off-axis angles. 



eter estimation procedure as the real data, are required. 
The simulations in H93 demonstrated that source confu- 
sion sets the ultimate limit in ROSAT deep survey work 
with the PSPC. The crowded-field multi-ML detection 
algorithm used in the current paper was specifically de- 
signed to better cope with source confusion. Therefore 
we felt it necessary to calibrate its efficiency and verify 
its accuracy through new simulations. We have simulated 
sets of PSPC and HRI observations with 200 ksec and 
lMsec exposure, respectively, approximating our current 
observation times. In order to compare our results with 
those of shallower PSPC surveys in the literature we sim- 
ulated additional sets of observations of 50 and 110 ksec 
each. The soft X-ray log(N)-log(S) function (0.5-2 keV) 
has been simulated according to the ROSAT deep sur- 
vey findings (H93) as a broken power law function for the 
differential source counts (all fluxes are given in units of 
10~ 14 erg cm~ 2 s _1 ): 

n(S) = ni x S~ bl for S > S b , n x = 238.1 b x = 2.72 

= n 2 x S~ b2 for S < S b , n 2 = 111.0 b 2 = 1.94 (3) 
S b = 2.66 

As shown in H93, the final results are relatively indepen- 
dent from the actual log(N)-log(S) parameters chosen for 
the artificial fields. In the simulations point sources are 
placed at random within the FOV, with fluxes drawn at 
random from the log(N)-log(S) function down to a mini- 
mum source flux of 4.65 x 10~ 18 , where formally all the X- 



ray background is resolved for the assumed source counts. 
(For the later comparison between input and output cat- 
alogues only input sources with fluxes larger than 10~ 16 
or more than 5 simulated photons are maintained in the 
input catalogue.) For each source the ROSAT vignetting 
correction, the corresponding ECF (see Table 2) and the 
exposure time (50, 110, 200 ksec for PSPC, lMsec for 
HRI) are applied to the source flux to obtain the expected 
number of source counts. The actual source counts are 
drawn from a Poissonian distribution and folded through 
the point spread function. The realistic multi-component 
point spread function model is taken from Hasinger et al. 
(1994b) for the PSPC and from David et al. (1996) for the 
HRI. Finally, all events missing in the field, i.e. particle 
background, non-source diffuse background and photons 
not simulated in the log(N)-log(S) function are added as 
a smooth distribution to the image. 

A total of 66 PSPC fields of an exposure of 200 ksec, 
50 fields of 50 ksec and 27 fields of 110 ksec each has 
been simulated in the 0.5-2 keV band. The corresponding 
number of HRI fields with 1 Msec exposure was 100. All 
simulated exposures have been analysed through exactly 
the same sequence of detection, background estimation 
and parameter estimation algorithms as the real data. The 
200 ksec PSPC images were analysed in two ranges of off- 
axis angles: 0-12.5 arcmin and 12.5-18.5 arcmin, the HRI 
lMsec images for off axis angle ranges of 0-10 arcmin and 
10-15 arcmin and the 50 and 110 ksec PSPC exposures 
inside 15 arcmin. 
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6.1. Flux conservation 

For each detected source the process of "source identi- 
fication" has been approximated by a simple positional 
coincidence check. A detected source was identified with 
the counterpart from the input list, which appeared clos- 
est to the X-ray position within a radius of 30 arcsec. 
A direct comparison between the detected and simulated 
number of photons for sources with fluxes brighter than 
5 x 10~ 14 cgs shows, that the scale of the output flux is 
the same as that of the input flux to an accuracy of better 
than 1%. This means that at bright fluxes no systematic 
errors exist, apart from the PSF loss factor of 0.90 (see Ta- 
ble 2), which is corrected for globally in the simulations 
as well as in the real data. This correction is necessary, 
probably due to the misfit between the broader compo- 
nents of the real PSF and the simple Gaussian model in 
the multi-ML fit. For the HRI this PSF loss factor is 0.92 
(see Table 2). 

6.2. Limiting sensitivity and logfN)-log(S) relation 

A log(N)-log(S) function has been derived for the simu- 
lated detected sources (in rings of 0-f 2.5 arcmin and 12.5- 
18.5 arcmin for the PSPC and rings of 0-10 arcmin and 
10-15 arcmin, respectively for the HRI). For easier com- 
parison the differential source counts n(S) were divided by 
the source counts expected for a Euclidean distribution, 
i.e.: 

n{S) cor = n{S) x S 2 - 5 (4) 

Fig. H shows the comparison between the input log(N)- 
log(S) and the detected number counts in this representa- 
tion for both 200 ksec PSPC and the 1 Msec HRI simula- 
tions. The faintest sources detected in the PSPC at small 
off-axis angles have a flux of 2 x 10~ 15 cgs. For a flux lower 
than 3 x 10 -15 cgs, the source counts fall significantly be- 
low the simulated log(N)-log(S) function, while above this 
flux there is a slight excess. These effects are well known 
from ROSAT PSPC deep survey simulations and are due 
to a number of selection and confusion effects. At larger 
off-axis angles the sensitivity is reduced correspondingly. 
The faintest HRI sources reach down to fluxes of 10 -15 cgs 
and the deviation between input and output source counts 
becomes insignificant above a flux of 1.8 x 10~ 15 cgs for 
small off-axis angles. Again, at larger angles the sensitiv- 
ity is reduced. The source confusion effects are much less 
pronounced in the HRI. 

The deviation between input- and output source 
counts in the simulations can be used to correct the ob- 
served log(N)-log(S) function down to the faintest limit- 
ing fluxes (see H93). Utilizing all sources detected in the 
real PSPC and HRI pointing observations inside off-axis 
angles of 18.5 and 15 arcmin, respectively, we calculated 
the corrected observed log(N)-log(S) function displayed in 
Fig. HI The new PSPC data agree very well with the source 



counts published earlier and extend those down to a flux 
of 2 x 10 -15 cgs and a surface density of 640 ± 75 deg~ 2 . 
The HRI log(N)-log(S) function for the first time extends 
the source counts down to fluxes of 10 -15 cgs and reaches 
a surface density of 970 ± 150 deg~ 2 , about a factor of two 
higher than any previous X-ray determination of source 
counts. All data are consistent with the previous determi- 
nation of source counts and fluctuation limits (H93). 



ROSAT logN-logS function 
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Fig. 3. Measured ROSAT log(N)-log(S) function in the Lock- 
man Hole. Filled circles give the source counts from the 207 
ksec PSPC observation. Open circles are from the ultradeep 
HRI observation (1112 ksec), slightly shifted in flux in order to 
avoid overlap of error bars. The data are plotted on top of the 
source counts (solid line) and fluctuation limits (dotted area) 
from Hasinger et al. (1993). The dotted line at bright fluxes 
refers to the total source counts in the RIXOS survey (Mason 
et al., 1996, priv. comm.). 



6.3. Source confusion 

The detected source catalogues are affected by biases and 
selection effects present in the source-detection procedure. 
The most famous of those is the Eddington bias, which 
produces a net gain of the number of sources detected 
above a given flux limit as a consequence of statistical er- 
rors in the measured flux (see discussion in H93). Another 
selection effect, most important in the deep fields consid- 
ered here, is source confusion. The net effect of source 
confusion is difficult to quantify analytically, because it 
can affect the derived source catalogue in different ways: 
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Fig. 4. Detected flux versus input flux for (a) 66 simulated PSPC fields of 200 ksec exposure, and (b) 100 simulated HRI fields 
of 1 Msec exposure each. The PSPC data is for off-axis angles smaller than 12.5 arcmin, for the HRI the limit is 10 arcmin. 
Fluxes are in units of 10 -14 erg cm~ 2 s _1 . 



(1) two sub-threshold sources could be present in the 
same resolution element and thus mimic a single detected 
source. This leads to a net gain in the number of sources, 
similar to the Eddington bias. 

(2) two sources above the threshold could merge into a 
single brighter source. In this case one source is lost and 
one is detected at a higher flux. Whether the total flux 
is conserved or not depends on the distance between the 
two sources and on the details of the source detection al- 
gorithm. 

(3) the detection algorithm cannot discriminate close 
sources with very different brightness, which results in a 
net loss of fainter sources. 

The effects of confusion become immediately obvious 
in fig [| where for each detected source the detected flux 
is compared to the flux of the nearest input source within 
30". While for bright sources there is an almost per- 
fect match, there is a significant systematic deviation for 
fainter X-ray sources, where most objects appear at fluxes 
significantly brighter than their input counterparts, which 
is a direct indication of confusion because every detected 
source is only associated with one input source while its 
flux may be contributed from several sources. 

We can quantify "confusion" e.g. by looking for 
sources, whose detected X-ray flux is significantly larger 
than their input flux. This means that a substantial frac- 
tion of the detected photons originates from other, con- 
taminating sources. We define SR as the ratio of the de- 



tected flux and the input flux increased by 3 times the 
statistical error: 



Sinp l+3fJS 

where as is the fractional flux error. For practical pur- 
poses a flux ratio of SR > 1.5 has been chosen, above 
which significant position deviations might be expected. 
In Table 5, we give the fraction of contaminated sources 
FC with flux ratio SR larger than 1.5 as a function of 
detected flux for different surveys. Significant confusion 
(larger than 10%) sets in below a flux of 5 x 10 -15 cgs for 
the 200 ksec PSPC observation and below 10 -14 cgs for 
the 50 ksec observation. In contrast, the HRI confusion 
is negligible over almost the complete range of detected 
fluxes. 

6.4- Effect of source confusion on optical identification 

While it is obviously possible to correct for confusion ef- 
fects on the source counts in a statistical way, the optical 
identification process relies on the position of individual 
sources. An independent measure of confusion is the frac- 
tion of sources with significant positional deviation from 
their respective counterparts. Fig. ^| shows for each sim- 
ulated source the distance to the nearest output source 
as a function of simulated input flux. As expected, bright 
sources have relatively small positional errors and there- 
fore small deviations of a few arcsec between simulated 
and output position. For fainter sources there is an in- 
completeness limit. More and more input sources are not 
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Table 5. Confusion estimates" 
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a Fraction of confused sources in percent according to different confusion measures (see text): FC is the percentage of contami- 
nated sources, having an output flux significantly larger than their input flux (SR > 1.5). FL is the percentage of lost sources, 
having no detected source within 15" from the input position. FU is the fraction of unidentifiable detected sources, having no 
input source within 15" from the detected position. b detected flux for FC and FU, input flux for FL 



PSPC 200 ksec 
Simulation 




HRI 1 Msec 
Simulation 




Fig. 5. Position deviation between input source and nearest output object as a function of input flux for the simulations in 
Fig. 0. The plot is cut off artificially for distances larger than 200 arcsec. 



detected at all and therefore are matched to the nearest 
(but wrong) detected X-ray source. For the HRI simula- 
tion this yields a well-defined cloud of spuriously matched 
sources below a flux of ~ 2 x 10~ 15 cgs. For the PSPC 
this incompleteness limit sets in below ~ 4 x 1CP 15 cgs 
and the spurious match cloud is not disjunct of the de- 
tected source correlation track. Confusion leads to a halo 
of PSPC sources with significantly larger positional devi- 
ation than the bulk of detected sources at the same flux. 
A practical limiting radius for the optical identification of 
faint ROSAT X-ray sources is 15", which is used by sev- 
eral authors (see e.g. Bower et al., 1996, Georgantopoulos 
et al., 1996). In Table 5 we give the fraction of lost sources 
FL(> 15") which have no detected source within a dis- 
tance of 15" from the input position. 



The quantity FL is a good tool to define a reliable flux 
limit for a survey, because it combines both the confusion 
and incompleteness effects. For practical applications we 
need, however, a quantity that can be applied to existing 
surveys, i.e. to detected sources. To do this, we invert the 
question that led to the quantity FL and ask: how many of 
the detected sources appear at a detected position further 
away from any "reasonable" input position than e.g. 15". 
For this purpose we exclude as reasonable input positions 
sources with input fluxes more than a factor 3 below the 
detected flux. This gives us an estimate of the fraction of 
unidentifiable sources in a survey, which we denote with 
the quantity FU. 

The quantities FC, FL and FU are compared in Table 
5 for different surveys and flux ranges. This table shows 
that, as expected, confusion effects increase dramatically 
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towards fainter fluxes. For fluxes above the incomplete- 
ness limit the different confusion measures, FC, FL and 
FU lead to similar estimates for the number of confused 
sources. Below this limit, FL increases dramatically, indi- 
cating severe incompleteness of the survey in addition to 
confusion. Using this table we have chosen the relatively 
conservative limits in flux and off-axis angles employed 
in the current survey, i.e. a flux limit of 5.5 x 10~ 15 for 
off-axis angles smaller than 12.5 arcmin and a flux limit 
of 1.1 x 1CP 14 for off-axis angles between 12.5 and 18.5 
arcmin. In the discussion section we compare these simu- 
lations to shallower PSPC surveys in the literature. 

7. Discussion 

With a limiting flux of ~ 10 -15 erg cm~ 2 s _1 , the X-ray 
survey in the Lockman Hole represents the deepest X-ray 
survey ever performed. The total observing time invested 
is quite comparable to that of other major astronomical 
projects, like e.g. the Hubble Deep Field (Williams et al., 
1996). Because of the expected confusion in the PSPC 
it was clear from the beginning, that HRI data would be 
necessary to augment the PSPC identification process. Be- 
cause of the smaller field-of-view of the HRI and because 
of its lower quantum efficiency, it was necessary to in- 
vest about a factor five more HRI time than PSPC time. 
Both the HRI raster scan, which provides excellent posi- 
tions for all brighter objects in the PSPC field, and the 
HRI ultradeep survey in the (slightly offset) central part 
of the field allow almost complete optical identifications of 
sources down to 5.5 x 10~ 15 erg cmT 2 s _1 (see paper II). 
Apart from this, the ultradeep HRI data provides a survey 
in its own right, which is not yet fully exploited. It will 
ultimately lead to reliable optical identifications a factor 
of 3-4 deeper than the current survey. On the other hand, 
the Lockman Hole data arc also of fundamental impor- 
tance for other observations with the ROSAT HRI since 
they provided for the first time a significant determination 
of the HRI scale factor correction which is important for 
all high-quality astrometry with the HRI. 

7.1. The ROSAT log(N)-log(S) function 

The new PSPC and HRI data shown in Fig. | carry the 
resolved source counts a factor of 2.5 deeper than the 
most sensitive previous determinations (H93, Branduardi- 
Raymont et al., 1994). It is very reassuring to see that the 
source counts are still consistent with the previous fluc- 
tuation analyses (H93, Barcons et al., 1994). There are, 
however, some problems at the bright end of the source 
counts. In our previous paper we used the source counts 
from the EMSS (Gioia et al., 1990) to constrain the bright 

end (V 5 -2.ofc e y) > L8 x 10 ~ 13 ) of the log(N)-log(S) re- 
lation. The first study to note a discrepancy between the 
EMSS and ROSAT counts was the RIXOS survey. Page et 
al. (1996) found, that the EMSS log(N)-log(S) for AGN is 



significantly lower than the one from RIXOS. The source 
counts from the ROSAT All-Sky-Survey Bright Source 
Catalogue (Voges et al., 1996) confirm this trend and ex- 
tend it to fluxes well above 10~ 12 erg cm~ 2 s _1 (Hasinger 
et al., 1997, in prep.). It appears that the log(N)-log(S) 
function at fluxes brighter than 10 -13 erg cm~ 2 s _1 has 
a slope of Pq ~ 2.3 and a normalization of ~ 91 and 
is therefore significantly flatter than what we assumed in 
H93 from a fit to the EMSS total source counts. The flatter 
log(N)-log(S) slope at bright fluxes has some consequences 
for the calculation of the resolved fraction. 



7.2. The resolved fraction 

As in H93 we restrict the analysis of the resolved fraction 
to the energy band 1-2 keV in order to minimize galac- 
tic contamination in the X-ray background spectrum. The 
absolute level and the spectrum of the X-ray background 
in this energy range is still a matter of debate (see the dis- 
cussion in Hasinger 1996). In H93 we assumed a 1-2 keV 
background flux of 1.25 x 10 -8 erg cmT 2 s _1 sr _1 . From 
ASCA data in the energy range 1-7 keV, Gendreau et al. 
(1995) have derived a flat spectrum with photon index 
1.4 and a normalization of 8.9 keV cm~ 2 s _1 keV" 1 at 1 
keV. Integrating this spectrum we derive a lower limit to 
the 1-2 keV XRB flux of 1.22 x 1(T 8 erg cmT 2 s~ x sr _1 . 
This can be contrasted to the earlier ROSAT PSPC de- 
termination of a steep spectrum in the 0.5-2 keV band, 
with photon index around 2 and a normalization around 
13 (see Hasinger 1996). Integrating this over the 1-2 keV 
band yields a flux of 1.44x 10 -8 erg cmr 2 s _1 sr _1 , which 
is consistent with the background flux derived by Chen et 
al., 1997 from a joint ROSAT/ASCA fit (1.46 in the same 
units). 

In H93 we integrated the analytic log(N)-log(S) 
function (Eq. 3) for fluxes brighter than 2.5 x 
10 -15 erg cm~ 2 s _1 , including a fit to the EMSS at fluxes 
brighter than 1.8 x 10~ 13 erg cmT 2 s _1 . We obtained a 
resolved 1-2 keV flux of 0.74 x 10~ 8 erg cmT 2 s _1 sr -1 . 
If we follow the same prescription, but now integrat- 
ing the H93 source counts above our limiting flux of 
10 -15 erg cmT 2 s _1 , we arrive at a resolved flux of 0.89 in 
the same units. If we, however, correct for the flatter slope 
of the bright end log(N)-log(S) (see above), we arrive at 
a resolved flux value of 0.99 (same units). We use this re- 
solved flux in comparing with the total diffuse background 
in the 1-2 keV band. 

If we adopt a 1-2 keV background flux of 1.45 x 
10 -8 erg cm~ 2 s _1 sr~ 1 , we have resolved 68% of the 1-2 
keV X-ray background at a flux of 10~ 15 erg cm~ 2 s^ 1 . 
If, however, the lower ASCA spectrum with a flux of 
1.22 x 10 -8 erg cmT 2 s _1 sr -1 holds, we have already re- 
solved 81% of the background. We see, that the uncer- 
tainty in the resolved fraction is now dominated by the 
systematic error in the absolute background flux and not 
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Table 6. Comparison of ROSAT surveys 



Survey" Area b T c Sf im Off° Err f Sf. ange FC h FL 1 FU j 



RIXOS 


~ 20 


>8 


3.0 














CRSS 


3.9 


>6 


2.0 














DRS 


1.4 


21-49 


0.3 


15 


15 


0.3-0.6 


13 


62 


23 


MARA 


0.2 


55 


0.5 


15 


15 


0.5-1 


11 


33 


19 


NEP 


0.2 


79 


1.0 


15.5 


15 


1-2 


6 


6 


7 


UKDS 


0.2 


110 


0.2 


15 


10 


0.2-0.4 


10 


76 


42 


RDS 


0.3 


207 


0.5 


12.5 


15 


0.5-1 


7 


10 


5 


HRI 


0.1 


1000 


0.2 


12.5 


5 


0.2-0.4 


2 


7 


5 



a cf. text for explanation of acronyms 
b approximate, in square degrees 
c PSPC exposure time, in ksec 

d quoted limiting flux £(0.5 — 2kcV) in units of 10~ 14 erg cm~ 2 s _1 

8 maximum off- axis angle in survey in arcmin 

f radius of error circle searched in survey in arcsec 

g flux range for which FC, FL and FU have been estimated 

h percentage FC of contaminated objects 

1 percentage FL of lost objects 

J percentage FU of unidentifiable objects 



by the source counts. As a best guess for the resolved frac- 
tion we take 70 - 80%. 

7. 3. Comparison to shallower PSPC surveys 

The optical identifications in the Lockraan Hole and the 
detailed comparison with other work is described in paper 
II (Schmidt et al., 1997). Here we want to draw global com- 
parisons with other PSPC surveys, in particular applying 
the simulation results described above. Quite a number of 
groups are involved in optical identifications of deep and 
medium-deep survey fields observed with ROSAT (see Ta- 
ble 6) and already many papers have been written about 
the interpretation of ROSAT survey results. Among the 
most debated findings is the claimed detection of a possi- 
ble new class of X-ray active, optically relatively normal 
emission line galaxies at faint flux levels (Griffiths et al., 
1996, McHardy et al., 1997). Unfortunately, however, only 
very few surveys so far have been published formally (i.e. 
including a catalogue of detected sources and a detailed 
description of the detection and identification procedure), 
so that quantitative comparisons can be made. Among 
those published are the Cambridge-Cambridge ROSAT 
Serendipity Survey (CRSS, Boyle et al., 1995), a small 
part of the deep ROSAT survey (DRS, Shanks et al., 1991, 
see also Gcorgantopoulos et al., 1996) and the North Eclip- 
tic Pole Survey (NEP, Bower et al., 1996). Among those 
waiting to be fully published are the ROSAT International 
X-ray Optical Survey (RIXOS, Mason et al., in prep., see 
also Page et al., 1996), the UK ROSAT deep field survey 
(UKDS, McHardy et al., 1997), the full DRS, the survey in 
the Marano field (MARA, Zamorani et al., in prep.), and 



finally our own ROSAT Deep Survey (RDS, this paper; 
Schmidt et al., 1997). Nevertheless, the global properties 
of these surveys are known and some details can be ob- 
tained from the existing literature. See Table 6 for a sum- 
mary of the surveys in question. The table also includes a 
prediction for the HRI ultradeep survey. 

Using the experience gained from the ROSAT Deep 
Survey and applying the simulation tools to the shallower 
surveys we can now make some statements about the ex- 
pected quality of the other surveys, in particular with re- 
spect to possible confusion and the corresponding opti- 
cal misidentification. As we have seen in Table 5, vari- 
ous confusion problems, as indicated by the FL, FC and 
FU percentages, become severe roughly within a factor 
of 2 from the formal detection threshold of a PSPC ex- 
posure longer than 50 ksec. It is interesting to note that 
this behaviour is not a strong function of the exposure 
time. In order to obtain a more quantitative assessment 
of the different surveys, we use our PSPC simulations of 
50 ksec (DRS, MARA), 110 ksec (UKDS) and 200 ksec 
(RDS) exposures applying as far as possible the detailed 
information about exposure times, quoted flux limit, off- 
axis angles and assumed error circle radii in the individual 
surveys (see Table 6). For the 80 ksec NEP survey we use 
quantities interpolated between the 50 ksec and the 110 
ksec simulations. For lack of even shallower simulations 
we do not make statements about the RIXOS and CRSS 
surveys. For each of the deep and medium-deep surveys 
in Table 6 we calculate the quantities FC, FL and FU, 
for a flux range within a factor of two from the sensitivity 
limit quoted by the authors. This way we can predict for 
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each survey the fraction of misidentified sources among 
the faint source population, where presumably the new 
discoveries would be expected. 

It is not surprising to see, that those surveys that 
employ a relatively conservative limiting sensitivity (as 
judged from the ratio between exposure time and flux 
limit) have very small fractions of unidentifiable sources. 
This is confirmed by the high rate of identifications e.g 
for the NEP survey (Bower et al., 1996) and the RDS 
(Schmidt et al., 1997). The RDS is additionally helped 
by the HRI data in the field. Larger fractions of uniden- 
tifiable sources on the order of 10% are expected for the 
DRS and the Marano survey. Indeed, in the Marano field, 
which unfortunately does not have HRI coverage, we see 
a significant number of empty error boxes (Zamorani et 
al., in preparation). Unidentified source fractions of order 
20% are no major problem as long as one is studying ma- 
jority populations of X-ray sources (e.g. AGN). They are, 
however, already a substantial problem if one tries to iden- 
tify new classes of objects which necessarily are minority 
classes. There the error can easily approach 100%. For the 
UKDS, which tried to push deepest in the optical identi- 
fication, unfortunately the largest unidentifiable fraction 
(42%) is expected. Indeed, about 26% out of a total of 34 
UKDS sources fainter than 5 x 10~ 15 erg cm~ 2 s^ 1 are 
unidentified (as derived from Fig. 5 in McHardy et al., 
1997), either because the error boxes are empty or con- 
tain too faint objects. Judging from our simulations we 
also expect some of the UKDS sources with proposed op- 
tical counterparts to be misidentified. 

8. Conclusions 

We have presented the complete catalogue of 50 
X-ray sources with 0.5-2 keV fluxes above 5.5 x 
10~ 15 erg cm~ 2 s _1 from the ROSAT Deep Survey point- 
ing of exposure time 207 ksec in the Lockman Hole. The 
X-ray positions are largely defined by additional ROSAT 
HRI observations of more than lMsec total exposure time. 
Using the HRI data we have derived a new log(N)-log(S) 
function reaching a source density of 970 ± 150 deg~ 2 at 
a limiting flux of about 10~ 15 erg cmT 2 s" 1 . At this level 
70-80% of the 0.5-2 keV X-ray background is resolved into 
discrete sources. The observations and the analysis proce- 
dure specifically developed to cope with confused X-ray 
observations have greatly profited from detailed simula- 
tions of PSPC and HRI fields. Based on these simulations 
we have defined conservative limits on flux and on off-axis 
angles, which guarantee a high reliability of the catalogue. 
The fraction of confused or unidentifiable sources in the 
catalogue is expected at only a few percent, but is proba- 
bly even lower due to the existence of the deep HRI data. 

We also discussed simulations of shallower fields and 
show that surveys, which are only based on PSPC ex- 
posures larger than 50 ksec, become severely confusion 
limited typically a factor of 2 above their 4a detection 



threshold. Sizeable fractions of confused and unidentifiable 
sources are expected for some of the published surveys. 
This may have consequences for some recent claims of a 
possible new source population emerging at the faintest 
X-ray fluxes. 
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